Prediction of thickness limits of ideal ultrathin polar films 
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Competition between electronic and atomic reconstruction is a constantly recurring theme in 
transition metal oxides. We use density functional theory calculations to study this competition 
for a model system consisting of a thin film of the polar, infinite- layer structure ACuCb (A = Ca, 
Sr, Ba) grown on a non-polar, perovskite SrTi03 substrate. A transition from the bulk planar 
structure to a chain-type thin film accompanied by substantial changes to the electronic structure 
is predicted for a SrCuC^ film fewer than five unit cells thick. An analytical model explains why 
atomic reconstruction becomes more favourable than electronic reconstruction as the film becomes 
thinner and suggests that similar considerations should be valid for other polar films. 

PACS numbers: 68.35.-p, 68.35.Ct, 73.20.-r 



Introduction: The lowering of symmetry at surfaces 
and interfaces frequently results in new or enhanced phys- 
ical properties. A central goal of transition metal oxide 
thin film engineering is to exploit this by making interface 
effects dominate bulk properties in a controlled fashion 
[1, 2]. When a polar thin film is grown on a non-polar 
substrate, achieving this control can be very difficult be- 
cause an instability arises that can drive electronic and 
atomic reconstruction. In the case of LaA10 3 (LAO) 
thin films grown on SrTi0 3 (STO) substrates [2-5], the 
alternate stacking of positively (LaO + ) and negatively 
(A10^~) charged layers should result in huge internal elec- 
tric fields. If nothing else were to happen, the increas- 
ing electrostatic energy would give rise to a "polar elec- 
trostatic instability". The displacement in response to 
this instability [5, 6] of electrons (charge transfer or leak- 
age) and ions (atomic relaxation and reconstruction) can 
generate compensating electric fields. Whereas charge 
transfer and atomic relaxation result in atomically sharp 
interfaces and ideal thin films with essentially intrinsic 
properties [7-13], atomic reconstruction can result in the 
formation of structures with entirely different intrinsic 
properties, or in changes to the local stoichiometry due 
to e.g. oxygen vacancy formation [14] or ionic intermix- 
ing [15, 16]. The competition between electronic charge 
transfer and atomic relaxation and reconstruction in ox- 
ide thin films is the source of much puzzling behaviour 
[5, 14-16]. 

To unravel the details of this competition is a chal- 
lenge for experiments because of the limited resolution 
of interface sensitive measurements and the low visibil- 
ity of oxygen atoms in most techniques [5]; it is also a 
challenge for theoretical studies to describe atomic recon- 
struction that results from subtly different experimental 
conditions [13, 17] when the role of stoichiometry is un- 
clear. In this paper, we use first-principles total energy 
calculations to study the growth of thin films of the po- 
lar infinite- layer copper oxide ACu02 (A = Ca, Sr, Ba) 
on a non-polar perovskite STO substrate and predict a 



stoichiometric atomic reconstruction as a function of the 
film thickness. As the parent compound of cuprate high- 
temperature superconductors (HTS), ACu02 (ACO) has 
been intensively studied [18-21]. The recent discovery of 
very high mobilities in heterostructures containing thin 
films of SrCu02 (SCO) [22] and a theoretical proposal for 
electron-hole liquids [23] has led us to reexamine its struc- 
tural and electronic properties. Our finding of chain- 
type formation in thin films as a result of an electro- 
static instability and atomic reconstruction provides im- 
portant insight into the growth of thin films of cuprate 
HTS. We note that the chain- type formation is differ- 
ent from a GdFe03-type [9] or a polar [11] distortion in 
LAO | STO, ferroelectric or antiferrodistortive distortions 
in PbTi0 3 |SrTi0 3 [24] and in bulk Ca 3 Mn 2 7 [25]. 

Method: We study ACO thin films in three forms: 
grown on an STO substrate; in ACO | STO multilayers; 
and freestanding. The thickness of ACO is varied from 
one to six unit cells while keeping the STO thickness fixed 
at five unit cells. The in-plane lattice constants of all thin 
films are fixed at the equilibrium value of the bulk STO 
substrate calculated to be asTO = 3. 87 A and all atoms 
are allowed to relax fully. The local density approxima- 
tion (LDA) calculations were carried out with the projec- 
tor augmented wave method [26] as implemented in the 
Vienna Ab-initio Simulation Package (VASP) [27]. Us- 
ing the generalized gradient approximation (GGA) does 
not change our main conclusions. To better describe lo- 
calized Cu d electrons, we use the rotationally invariant 
LDA+U method [28] with U - J = 6.5 eV [29]. 

Bulk materials: Bulk STO is a well-studied band in- 
sulator and a popular substrate for growing thin oxide 
films. It has a typical perovskite structure in which Sr 
and Ti06 form a CsCl lattice and the Ti06 unit consists 
of an oxygen octahedron with Ti at its center and the 
oxygen atoms at the centers of the Sr cube faces. Assign- 
ing the formal ionic charges Sr 2+ , Ti 4+ and O 2- , it can 
also be described as an alternate stacking of uncharged 
SrO° andTiO^ layers [2]. 
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FIG. 1. Left panel: bulk SrCu02 with infinite- layer structure. 
Oxygen atoms are marked by blue spheres. Red and green 
represent surfaces of constant magnetization density but op- 
posite sign. Right panel: the red curve illustrates the strain 
energy of bulk SrCuCb on varying the lattice constant in the 
c direction keeping the in-plane lattice constants a and b fixed 
at asTO- When the structure is rotated so that b and c are 
forced to match asTO and a is varied, the resulting strain 
energy is given by the blue curve. 



The copper oxide AC11O2 with infinite layer structure 
[18] depicted in Fig. 1 can be regarded as a defect per- 
ovskite with ordered oxygen vacancies. It consists of pos- 
itively charged A 2+ and negatively charged CuO 2- lay- 
ers. Because all of the oxygen sites in the AO plane are 
vacant, there are Cu-0 bonds in the xy plane but none 
in the z direction. The missing atoms and bonds lead 
to a reduction of the lattice constant in the z direction, 
reducing the symmetry from cubic in STO to tetragonal 
in ACO. The lattice constant of ACO increases as the ion 
radius increases in the series A = Ca, Sr, Ba. The calcu- 
lated values given in Table I agree well with experiment 
and SCO is seen to have the smallest lattice mismatch 
with STO. Though a complex orthorhombic structure for 
SCO can also be synthesized [30], only the infinite layer 
structure is taken into consideration in thin film growth 
experiments [21-23, 31], thus we refer here to the infinite 
layer structure as the bulk structure and neglect other 
possibilities in this paper. 

Including an on-site Coulomb repulsion term U only 
slightly affects the calculated lattice constant, and can re- 
produce the observed antiferromagnetic insulating (AF- 
I) ground state [32] with an energy gap located be- 



TABLE I. Lattice constants of perovskite SrTiOs and infinite 
layer copper oxides ACu02(A=Ca, Sr, Ba). 





a LDA 


c LDA 


a GGA 


C GGA 


a exp 


c exp 


SrTiOs 


3.87 


3.87 


3.95 


3.95 


3.905 


3.905 


CaCu0 2 


3.77 


3.08 


3.87 


3.20 


3.853 a 


3.177 a 


SrCu0 2 


3.84 


3.38 


3.95 


3.47 


3.926 b 


3.432 b 


BaCu0 2 


3.92 


3.68 


4.03 


3.84 








a Ref.20 
b Ref.19 



FIG. 2. Relaxed structures and magnetic isosurfaces of 
SrCu02|SrTi03 multilayers with one unit cell of planar- (left 
panel) or chain-type (right) SrCu02 thin film. For simplic- 
ity, only three unit cells of SrTiOs are shown. The arrows in 
the left panel schematically indicate the oxygen displacments 
that transform the planar into a chain-type structure. 



tween filled oxygen p bands and an unfilled minority spin 
Cu d x 2_ y 2 band [29]. Because its nodes point toward 
neighbouring oxygen atoms, the d x i_ y i orbital hybridizes 
strongly with the oxygen p states and is pushed up in en- 
ergy. The resulting Cu 2+ d 9 configuration corresponds to 
all Cu d orbitals being filled except for the minority spin 
d x 2_ y 2 orbital. This is apparent in the surfaces of con- 
stant magnetization density plotted in Fig. 1 where the 
spin density on the oxygen atoms also illustrates the AF 
superexchange coupling between neighbouring Cu ions. 
The hybridization is mainly confined to the Cu02 plane 
and does not affect the ionic character of the A 2+ and 
CuO^" 2 layers so the electrostatic instability of ACUO2 
thin films still exists by analogy with LAO thin films. 

ACu0 2 thin films: When SCO is grown on an STO 
substrate, it is assumed that positively charged A 2+ and 
negatively charged CuO 2- layers alternate [21]. Because 
this planar ("CUO2") structure will lead to an electro- 
static instability, we also consider a chain type ( "CuO" ) 
structure that is formed by moving one oxygen atom 
from the CuO 2- layer to the oxygen vacancy position 
in the A 2+ layers as illustrated in Fig. 2 for the case of 
SCO I STO multilayers. This results in a thin film consist- 
ing of uncharged SrO° and chain-type CuO layers that 
does not suffer from an electrostatic instability but at the 
expense of strain (Fig. 1). For sufficiently thick films, we 
expect the planar structure to be most stable. It is not 
clear a priori what will occur for thin films of thickness 
d. 

To answer this question, we calculate the energy dif- 
ference, AE = £ chain - £P lanar as a function of d. As 
shown in Fig. 3, AE depends strongly on both d and 
on the cation A. For SrCu02 thin films grown on a lat- 
tice matched STO substrate, it increases gradually from 
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FIG. 3. Total energy difference between planar and chain- 
type ACO thin films as a function of the film thickness d. 
Symbols represent LDA results for SrCuCb thin films in three 
forms: grown on a SrTiOs substrate (black circles), as a mul- 
tilayer with SrTiOs (squares), and freestanding (triangles); 
the unfilled squares represent LDA+U results for multilayers. 
The solid black line for SrCuCb is based on the analytical 
model including the electrostatic energy and strain energy 
discussed in the text while the red and blue solid lines differ 
only in terms of the calculated strain energies of CaCuCb and 
BaCu0 2 . 



-1.2 to +0.05 eV per unit cell, as the thickness increases 
from one to six unit cells. Since both structures have 
the same stoichiometry, the sign of AE is a direct mea- 
sure of their relative stability. Below a critical thickness 
of about five unit cells, chain type thin films are ener- 
getically favourable; above it, planar type films are more 
stable. In the cation series Ca— ^Sr— ^Ba, AE(d) decreases 
but depends on thickness in the same way for all ACu02 
thin films. Similar behavior is found for multilayers and 
freestanding thin films so surface or interface effects are 
not essential for a qualitative understanding. We then 
decompose AE(d) as Ec — Ep, where Ep (E p ) denotes 
the electrostatic energy (density) induced by the insta- 
bility of planar type thin films and Ec is the chemical 
bonding energy caused by the oxygen atom moving to 
form a strained chain-type structure. 

The planar thin films are terminated with positively 
charged Sr 2+ and negatively charged CuO?, - surfaces 
with nominal areal charge densities ±<r = ±e/a 2 where 
a is the in-plane lattice constant. They can be mod- 
eled as parallel plate capacitors [5, 6, 12, 13, 17, 23, 33] 
in which the plates are separated by a thickness d of 
bulk material with dielectric constant s. The electro- 
static energy density associated with the electric field 
a/eeo between the plates is E® = a 2 /(2eeo) or approx- 
imately 20/e eV per unit cell, where eo is the dielec- 
tric constant of vacuum. This constitutes a strong in- 
stability that will drive charge transfer that depends on 
the thickness and electronic structure of the thin films 
[23] . If we assume a rigid flat band approximation with 
constant density of states a, the transferred charge is 
<j/[l + a 2 eeo/(ad)]. As a result of the screening by the 



transferred charge, the residual electrostatic energy den- 
sity becomes E p (d) = E®/[1 + ad/(a 2 eeo)] 2 , where the 
screening factor clearly depends on a and thin film thick- 
ness d. For a large gap insulator with a = 0, E p is a 
constant [38]. For a good metal with a large a, E p ap- 
proaches zero. For a small band gap semiconductor with 
a moderate a, E p vanishes when d is large, consistent 
with bulk limits; if d is reduced, less charge will be trans- 
ferred and consequently E p will be enhanced. When the 
films are very thin, charge transfer cannot take place to 
quench the instability. 

In that case, Ep can be quenched by displacing oxy- 
gen atoms to form a chain-type structure as indicated in 
Fig. 2. Because the resulting structure corresponds to a 
strained ACu02 thin film grown in a (100) orientation on 
a (001) STO substrate (rhs of Fig. 2), rather than with 
a conventional (001) orientation (lhs of Fig. 2), doing so 
will change the bonding energy Ec- The two types of 
thin film only differ in terms of strain. To estimate the 
strain energy of planar type thin films, we fix the in-plane 
lattice constant of bulk ACu02 at asTO and minimize the 
total energy with respect to c. For chain type thin films, 
we consider b and c as the in-plane lattice constants, fix 
them to be equal to asTO and minimize with respect to 
a, as shown in Fig.l. The energy difference per unit cell 
between the two minima is Ec = 0.45 eV for SrCu02, 
1.19 eV for CaCu0 2 , and -0.01 eV for BaCu0 2 . 

As shown in Fig. 3, AE(d) can be fit very well using 
these values for Ec and values of a = O.le/V and e = 8 in 
the modified parallel plate capacitor model for Ep. The 
very good fit indicates that the competition between the 
two structures can be represented in terms of electronic 
and atomic reconstructions and that charge transfer can 
occur to partly screen the electrostatic instability without 
altering the atomic structures. However, as the films be- 
come thinner, the effect of screening decreases. For fewer 
than about five unit cells of SCO on STO, the widely as- 
sumed planar type thin film cannot be stabilized because 
of the instability. In spite of a considerable cost in strain 
energy, the chain type thin film should become energeti- 
cally favorable because of the lowering of the electrostatic 
polarization energy Ep. 

Atomic and electronic structure: To discuss the atomic 
and electronic properties of the novel chain-type thin 
films, we focus on the multilayers shown in Fig. 2 with a 
single unit cell of SCO sandwiched between five unit cells 
of STO. The parallel behaviour of the multilayer and thin 
film curves (black spheres and squares) in Fig. 3 indicates 
that the essential physics is the same so we can dispense 
with discussing the effect of dangling bonds at the free 
surface. In the chain-type film, there are stiff Cu-0 bonds 
in the growth direction. Because of this, the multilayer 
containing chain- type SCO is almost 0.5 A longer in the z 
direction than that containing planar SCO. This is illus- 
trated in Fig. 4(a) where the spacing between layers con- 
taining Sr atoms is plotted. Measurement of this spacing 
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FIG. 4. (a) Spacings between Sr layers for SC0|ST0 multi- 
layers with planar and chain-type monolayers of SCO relative 
to asTO for bulk STO. The black squares and line represent 
ideal stacking without epitaxial strain, (b) The Ti-0 or Cu-0 
corrugation in the ^-direction for each layer in both types of 
SCO | STO multilayer. 



capacitor model including charge transfer. As the polar 
film is made thinner, less charge can be transferred, the 
residual electrostatic energy increases, and formation of 
non-polar chain-type structures becomes an energetically 
favorable way to resolve the instability without changing 
the stoichiometry. Similar considerations can be applied 
to LAO | STO. Although pure charge transfer can avoid 
the "polar catastrophe" [5, 11], it can not fully quench 
the electrostatic instability. A large residual electrostatic 
energy, especially in the low thickness LAO regime, will 
drive the formation of oxygen vacancies or other defects. 
This work was supported by "NanoNed" , a nanotech- 
nology programme of the Dutch Ministry of Economic 
Affairs. This work is also part of the research programs 
of "Stichting voor Fundamenteel Onderzoek der Materie" 
(FOM) and the use of supercomputer facilities was spon- 
sored by the "Stichting Nationale Computer Faciliteiten" 
(NCF), both financially supported by the "Nederlandse 
Organisatie voor Wetenschappelijk Onderzoek" (NWO). 



should be possible and would provide experimental evi- 
dence for chain type thin films. The presence of zigzag 
Ti-0 buckling in Fig. 2 and 4(b) in the planar-type films 
is a strong indicator of internal electric fields [11-13, 17] 
that should be observable in structural studies. [39] 

The ground state of the planar-type SCO film is an 
AF metal. Consistent with a previous LDA+U calcula- 
tion for YBa2Cu306 interfaces with STO [37], we find 
that the d x i_ y i orbital character of the Cu 2+ d 9 minor- 
ity spin hole is largely unchanged by the polar nature of 
thin films and the presence of apical oxygen atoms in the 
STO substrate that leads to Cu-O-Ti type bonding. In 
contrast, the ground state of the chain- type film is an AF 
insulator with a localized magnetic moment of 0.6/i# on 
Cu ions and an energy gap of 1.6 eV. The energy gap 
between filled O p states and an unfilled minority spin 
Cu d state will decrease on increasing the thickness of 
the SCO film. The unfilled Cu d orbital character will be 
d z 2_ y 2 , as shown in Fig. 2, because the Cu ions are still in 
a CUO4 square planar configuration surrounded by four 
oxygen atoms in the y-z plane rather than the original, 
bulk x-y plane. Our results suggest an explanation for 
the d 3z 2_ r 2 orbital observed at the YBa2Cus07|LaMn03 
interface [7]; the electrostatic instability of Cu02 planes 
[23, 31] may be driving oxygen atoms (or vacancies) to 
move and effectively tilt the CUO4 square towards the 
z direction, favoring orbitals pointing in that direction. 
The polar nature of Cu02 planes may play a key role in 
the structure and properties of ultrathin films and inter- 
faces of HTS cuprates. 

Conclusion: Using first-principles calculations, we pre- 
dict that a polar electrostatic instability will lead to chain 
type structures being preferred when films of ACUO2 are 
sufficiently thin and interpret this with a parallel plate 
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